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Abstract 

We present a new result on the K + — ► 7r + 7r°7 decay measurement using stopped 
kaons. The best fit to the decay spectrum comprised of 10k events gives a branching 
ratio for the direct photon emission of [3.8±0.8{stat)±0.7(syst)] x 10" 6 in the vr+ ki- 
netic energy region of 55 to 90 MeV. This result has been obtained with the assump- 
tion that there is no component due to interference with the inner bremsstrahlung. 
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1 Introduction 



The chiral anomaly [1,2] is a fundamental property of chiral quantum field theories such 
as the Standard Model of quantum chromodynamics (QCD). Its theoretical origin and 
mathematical properties are well understood, while experimental tests, which are cru- 
cial for the theoretical basis of particle physics, are relatively rare. Although the chiral 
anomaly can be interpreted as a short-distant effect, it manifests itself most directly in the 
low-energy interactions of the pseudoscalar mesons [3-6]. The chiral anomaly manifests 
itself also in the non-leptonic weak interactions. As already shown in Refs. [7,8], only 
radiative K decays (among them K + — > 7r + 7r°7) are sensitive to this anomaly. 

The total amplitude of the K + — > 7r + 7r°7 decay can be written as a sum of inner 
bremsstrahlung (IB) and direct emission (DE) amplitudes. The IB amplitude is associated 
with the K + — > 7i + n° (K n2 ) decay, in which the photon is emitted from the outgoing 
7r + , and is completely predicted by the quantum electrodynamics in terms of the K n2 
amplitude [9] . The DE amplitude, in which the photon is emitted directly from one of the 
intermediate states of the decay, can in turn be decomposed into electric and magnetic 
parts. 

Within the framework of Chiral Perturbation Theory (ChPT), the magnetic amplitude 
appears first at 0{p A ). There are two different ways for its manifestation: the reducible 
amplitude [7], which can be derived directly from the Wess-Zumino-Witten functional [6], 
and the direct contributions [8,10], which are subject to some theoretical uncertainties. 
In both cases the direct emission amplitudes are dominated by the anomaly. For the 
electric amplitude there is no definite prediction from ChPT. The experimental study of 
the electric amplitude, which could be observed through an interference (INT) with the 
IB amplitude, is of high interest not only for ChPT but also for possible non-standard- 
model effects such as the CP-violation asymmetry between the K + — > 7r + 7r°7 and the 
K~ — > 7r~7r°7 decay widths. 

The remarkable feature of the K + — > 7r + 7r°7 and the Kl — > 7r + 7r~7, is that the 
normally dominant IB amplitude is strongly suppressed because it violates the A/=l/2 
rule. Therefore the DE amplitude is relatively enhanced. The DE component can also be 
isolated kinematically. These features simplify the experimental verification of the anoma- 
lous amplitude, making the K + — > 7r + 7T°7 decay an important instrument to investigate 
the low energy structure of QCD. 

A conventional way to express the differential decay rate of the K + — > 7r + 7r°7 decay is 
to use terms of T + and W , where T + is the kinetic energy of the ir + in the K + rest frame, 
and W is an invariant parameter defined as W = (P ■ q)(p+ ■ q)/(m^ + M^ + ). Here, P, p + 
and q are the 4-momenta of the K + , 7r + and 7, respectively. The differential decay rate 
can be written in terms of the IB component rate as [11] 

where A is the decay amplitude for the K + — > 7r + 7r°, and E and M are the electric and 
magnetic amplitudes of the DE component, respectively. 
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Table 1: Summary of the previous measurements of the direct photon emission in 
7r ± 7r°7 decays. 





Sort of 




Number 




Experiment 


experiment 


Kaon 


of events 


Br(DE) 


BNL 1972 [12] 


in-flight 


K± 


2100 


[1.56 ±0.35 ±0.5] x lO" 5 


CERN 1976 [13] 


in-flight 


K± 


2461 


[2.3 ±3.2] x 10" 5 


Protvino 1987 [14] 


in-flight 


K~ 


140 


[2.05 ±0.46i£!g] x 10" 5 


BNL E787 2000 [15] 


stopped 


K+ 


2 x 10 4 


[4.7 ±0.8 ±0.3] x lO" 6 


KEK E470 2003 [16] 


stopped 


K+ 


4434 


[3.2 ± 1.3 ± 1.0] x 10 -B 


Protvino 2004 [17] 


in-flight 


K~ 


930 


[3.7 ± 3.9 ± 1.0] x lO" 6 



The experimental status of the measurement of the DE component of the K + — > 7r + 7r°7 
decay is shown in Tabled] As one can see, there is a large discrepancy between the results 
of the first three [12-14] and the second three [15-17] experiments. These results can also 
be compared with corresponding theoretical predictions for the branching ratio of the DE 
component in the same kinetic energy region of the 7r + [7,8,10]. The above-mentioned 
situation encouraged us to perform a new analysis of the experimental data of the E470 
experiment. 

2 Experiment 

Experiment E470 was performed at the KEK 12 GeV proton synchrotron using the E246 
experimental apparatus constructed to search for T- violating transverse muon polarization 
in K + — > n fi + h / (K^ 3 ) decay [18]. In addition to the T- violation search, spectroscopic 
studies for various decay channels have also been successfully performed using this detector 
system [19]. The apparatus is shown in Fig. [TJ and is described in detail elsewhere [20-23]. 
A separated K + beam (tc + /K + ~ 7) of 660 MeV/c was used with a typical intensity of 
1.5 x 10 5 kaons per spill of 2 s duration with a 4 s repetition time. To distinguish K + 's 
from 7r + 's a Cherenkov counter [20] was used. A beam counter B0 [23], an assembly of 
22 plastic scintillating counters, placed before the Cherenkov counter, was used to obtain 
the K + and 7r + beam profiles. Kaons were slowed down in an Al±BeO degrader, and 
stopped in an active target made of 256 scintillating fibers, which was located at the 
centre of a 12-sector superconducting toroidal spectrometer. The charged ir + from kaon 
decay passed through one of the twelve fiducial (also called TOF1) counters, surrounding 
the target, and entered one of the spectrometer sector gaps, and was then detected by 
a TOF2 counter, located at the exit of the spectrometer. Tracking of charged particles 
and their momentum analysis was done using multi-wire proportional chambers at the 
entrance (C2) and exit (C3 and C4) of each magnet sector, as well as by the active target 
and an array of ring counters [21] surrounding the target. The energies and angles of 
the photons were measured by the CsI(Tl) photon calorimeter consisting of 768 CsI(Tl) 
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Figure 1: The layout of the detector for E470 experiment. 

modules [22]. The CsI(Tl) barrel had twelve holes for the charged particle entry to the 
toroidal spectrometer, and covered a solid angle of about 0.75 x4n sr. 
To collect the K + — > 7r + 7r°7 data the following trigger was used 

C K x [Fidi x TOF2i] x 37 x (1,2)50 x (l,2)ring, 

where Ck is a kaon Cherenkov counter, Fidi is a hit in the ith fiducial (TOF1) counter 
or in its immediate neighbors, TOF2i is a hit in the ith TOF2 counter, 37 is 3-photon 
clusters in CsI(Tl) barrel, and (1, 2)50 and (1, 2)ring are multiplicities of the BO counter 
and ring counters <2, respectively. 
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3 Analysis 



In this analysis we aimed to improve our previous result [16] of the measurement of the 
direct photon emission in K + — > 7r + 7r°7 decay by reanalyzing the experimental data ac- 
cumulated in 2001. For this purpose a new analysis strategy was adopted. This included: 
1) selecting as many good K + — > 7r + 7r°7 events as possible by imposing several loose 
cuts; 2) increasing the quality of data by improving the pairing efficiency of the photons 
from 7T° decay and 3) improving the K n3 background rejection. The analysis strategy 
was carefully studied and optimized by a Monte Carlo simulation. 

The decay time of the stopped K + in the target was defined by a signal in the TOF1 
counter. To reject in-flight K + decays this time was required to be more than 3 ns later 
than the K + arrival time into the target. 

For the analysis only events with three photon clusters in the CsI(Tl) calorimeter have 
been selected. This condition considerably suppresses the backgrounds from the K w2 , 
K e3 , K^ 3 and K n3 modes, although some fraction of these events pass the trigger due to 
different reasons- the K e3 and modes due to an accidental photon hit in the CsI(Tl) 
calorimeter; the K n3 mode due to the escape of one photon into a CsI(Tl) calorimeter 
hole. The K n2 events pass through the trigger condition mostly due to identification one 
of the two photons from n° decay dissipated in two parts as two photons. As a result of 
such misidentification we have three photon clusters in the CsI(Tl) calorimeter instead of 
two. 

To separate e + 's and /z + 's from 7r + 's, the squared mass of the charged particles M^ OF 
was calculated using the time of flight of the charged particles from the TOF1 to the 
TOF2 counter, and their specific energy deposition dE/dx in the TOF2. Such separation 
completely remove the remaining K e3 background, and considerably reduced the K^ 3 
background. 

Constraints on the reconstructed K + mass = E w + + Y%=i E-yi were set to be 

420 < Mk < 510 MeV/c 2 . To set a cut on the momentum px = p n + + Y^=\Pn we used 
the parameter Ap = ^Jp x 2 + p y 2 +p z 2 , where p a = p n+a + Ef=iP 7l a (« = x,y,z), and 
used the condition of Ap < 60 MeV/c. 

To reject in-flight decays of 7r + 's and suppress the scattering of charged particles 
from the magnetic pole faces, a track consistency cut \ 2 < 10 nas been imposed. The 
remaining K^i background has been removed by requiring the 7r + momentum to be less 
than 180 MeV/c. 

For further analysis we need to reconstruct the n°. Since we have three photons in 
the CsI(Tl), there are also three possible combinations to form the n°. To improve the 
pairing efficiency of the photons the following method was developed. The quantities Qf ■ 
(i, j=l,2,3), sums of six asymmetric terms for each pair of photons comprised of the ith 
and jth photon, were introduced as follows: 
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Table 2: The values of parameters in the Qf „• expressions. 



Expression 


Offset 


a (Expression <0) 


cr (Expression >0) 


E caic _ E ex P _ AE MeV 

76-i-j 76-i-j T' 

m 7i7j . - M 77 , MeV/c 2 
^o 7 . - 0,h> 7 , deg 
^f 7 % )e -^ 7j) -A^o,MeV 
6 ex P — Q-k+'y, deg 
0f+lo - ^+^o, deg 

»,.7 


A£ 7 = 2.54 
M 77 = 130.72 
^o 7 = 107.53 
AE n o = -2.58 
^+ 7 = 129.18 
^ +7r0 = 150.04 


= 6.25 
a mii = 22.05 
a 9o =19.52 

= 12.12 
a e , = 32.60 

7f-f 7 

a e + = 24.43 


a El = 12.86 

0"m 77 =4.5 
(T0 O =17.52 

CT^ = 9.31 

a e I = 15.10 

7r~r 7 

cr = 9.67 

7J--I- tj-U 



m. 



7i7j 



-M 77 ) 2 , (E^_ 3 



+ 



. -A£ 7 ) 2 (CV- W) ! 



o- £7 



(ere 



(2) 



Here, the superscripts "exp" and "calc" stand for the experimental and calculated values, 
respectively. So, Effi^.^ and E^ l ^. y .^ are the calculated energies of the 7r° comprised of 
a combination of zth and jth photons, using the formulae: 



rpcalcO 
*V>( 71 7i) 



\r2m w o 



\ 



v 1 COS ^7i,7J y 



1 - 



+ E exp 



(3) 



v 7j 



rricaicl 
£/ t°(7 ! 7 j ) 



M 2 , + m 2 n — m 2 , 2E?* (E ex + P - PT cos 6 ex + p 



2M 



2M, 



(4) 



calc 

76- 



the energy of the (6 



E, 

calculated from the following formula: 



j)th photon (supposing it is the unpaired one), is 



exp \ 



E: 



calc 
76-i- 



(M 2 K + m 2 n+ -m 2 n0 -2M K E e n 



(5) 



rn 



2E exv E exv (\ — cosOy^) is the invariant mass of the (i, j) photon combination. 
Other parameters in the Q\ • expressions, which were obtained using MC simulation, are 
given in Table EJ 

Since the distribution of the expressions in parentheses of each term's numerator in 
(2) is asymmetrical, the root mean square (a) corresponding to each term has two values, 
one corresponding to the expression being < 0, and the other to the expression being > 0. 
At this stage the variables of Qij have been optimized to maximize the pairing efficiency 
of the DE component using the MC simulation, which we obtained to be 88%. For the IB 
component we obtained 70% efficiency. Sub-indexes of the minimum Q it j are accepted to 
be the indexes of the two photons, which come from the ir°. We number all three photons 



6 



in the following way: the most energetic photon from the 7r° decay is ascribed the number 
1, the lower energetic photon is labeled by 2. We ascribe number 3 to the third unpaired 
photon. 

At the second stage we represent the events in a scatterplot as shown in Fig. 
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Figure 2: Scatterplot of 9 n +^ 2 vs E±f c0 . Here, 9 n + 72 is the opening angle between n + and 
72, and Eff c0 is the parameter Effi^^.^ with i = 1, j = 3. The detached region is a region 
of incorrect photons pairing. 



As the MC simulation shows the region in the rectangular box is a region of photons 
mispairing. Namely, for these events the 72 doesn't really come from the 7r° decay but is 
the unpaired one. Consequently the 73 is actually the lower energy photon from the ir° 
decay. In order to determine the true combination of the photons in the detached region as 
well, we renumber the events in this region as follows: we ascribe number two to the third 
photon and number three to the second one. Such simple cross-renumbering of events 
allowed us to enhance the pairing efficiency of the IB component up to 93%, although the 
DE component was reduced down to 81%. Nevertheless, on the whole we increased the 
pairing efficiency considerably, and this in turn reduced the systematic error associated 
with mispairing. For comparison, in the previous analysis [16] the pairing efficiencies were 
85% and 79% for IB and DE components, respectively. 

The K n3 background was rejected by a 2-dimensional cut as shown in Fig. El Such a 
cut not only allowed us to successfully reject this mode but also saved good K + — > n + n 0r y 
events in the 7r + low momentum region. 

The total number of good K + — > 7r + 7r°7 events in the tt + momentum region of 115 to 
180 MeV/c, obtained after applying the above selection criteria, is 10,154, which is about 
2.3 times higher than in the previous analysis [16]. The background from K n3 decay 
has been estimated to be less than 0.2% by the simulation. In the previous analysis the 
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Figure 3: Scatterplot of the 7r + momentum versus the sum of "invariant masses" comprised 
of 7i, 73 and of 72, 73. The shown line is a cut condition to remove the K n3 . 



contribution of this mode was about 1.2%. 

A MC simulation of the decay was used to test the simulation of the experiment. 
The matrix element of K n3 decay for this simulation was taken from the PDG [24]. In the 
simulation, the form factors of the K^a matrix element were taken to be the world average 
quoted by the PDG [24]. In this analysis we selected events with 4 photon clusters in 
the Cs(Tl) using the following conditions, < 135 MeV/c, 14000 < (M^°/) 2 < 25000 
(MeV/c 2 ) 2 , Ap < 80 MeV/c, 380 < M K < 520 MeV/c 2 , 80 < E n o < 200 MeV, 60 < 
E n o < 190 MeV, 60 < m {ll)1 < 150 MeV/c 2 , 55 < m (77)2 < 150 MeV/c 2 . Various K n3 
spectra are shown in Figs. 0J 03 and we see excellent agreement between the experimental 
data and the MC simulation. 

A MC simulation of the K^-y decay was carried out for both the IB and the DE 
components of the decay, using the matrix elements given in [25]. The detector 

acceptances for the IB and the DE components have been obtained to be Q(IB) = 
0.566 x 10 -3 and Vl(DE) = 1.659 x 10~ 3 . Figure M shows the K + — > 7r + 7r°7 spectra, which 
have been obtained after imposing an additional cut of W < 0.5 in order to remove the 
DE component. Again we find very good agreement between the experiment and the MC 
simulation. 
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Figure 4: K^s spectra: a) - 7r + momentum before correction for energy losses in the 
target; b) - tt + momentum after correction for energy losses in the target, c) - opening 
angle between n + and n®; d) - opening angle between n + and it®. The solid lines are the 
experimental data and the dashed lines are the MC predictions. 

4 Results 

The fraction of the DE component in the selected K + — > 7r + 7r°7 events has been obtained 
by fitting the experimental data distribution of O^+^o, E 1 and W in 3-D space with a sum 
of IB and DE components obtained from the MC simulation. The Maximum log-likelihood 
method has been used for fitting the \ 2 defined as follows 

X 2 = 2 E^ _n »+ n i ln ( n i/l/i)]. (6) 

where y^ — A x (y- B + ayf E ). Here, n« is the number of real data in the i-th bin of the 
histogram, and y\ B yf E are numbers of MC data of IB and DE components in the same 
i-th bin of the histogram, respectively. The variable a is the fraction of the DE component 
normalized to the IB one and A is a normalizing coefficient. The MINUIT package was 
used to obtain the minimum of \ 2 ■ The error values correspond to shifts of the variables, 
which increased the x 2 by 1. The value of a was obtained to be 0.026 ± 0.006 with a 
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Figure 5: Spectra of K n3 : a) - energy of E 7 x] b) - energy of E l2 ', c) - energy of E l3 ; d) - 
energy of E^. The solid lines are the experimental data and the dashed lines are the MC 
predictions. 

X 2 /ndf— 1.15. The experimental spectrum of W normalized to the pure IB component is 
shown in Fig. [7| The obtained branching ratio of the DE component in the 7r + kinetic 
energy region of 55 to 90 MeV is Br(DE) = [3.8 ± 0.8(stat)} x 10~ 6 . 

The major systematic errors come from the uncertainty of the CsI(Tl) calorimeter 
calibration, 7 mispairing effect, accidental background, the spectrometer field uncertainty, 
and also from the uncertainty of the kaon stopping x, y position. 

The uncertainty due to the 7 calibration of the CsI(Tl) calorimeter was estimated by 
changing the value of the DE branching ratio by the maximum shift of the calibration 
coefficients of the CsI(Tl) calorimeter. The corresponding shift of the DE branching ratio 
obtained in this way was ±0.5 x 10~ 6 . The 7 mispairing effect was studied by changing 
the pairing efficiency of the photons along with a change of the selection cuts. The shift 
of the DE branching ratio as a result of such changes was ±0.4 x 10~ 6 . The uncertainty 
from the accidental background, mainly associated with the decay, was estimated by 
varying the TDC gate widths, and was found to be less than 7% of the DE branching 
ratio value, that is ±0.3 x 10~ 6 . The error due to the spectrometer field uncertainty 
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o (tt + ,O (degree) E(j 3 ) (MeV) 

Figure 6: Spectra of K + — > 7r + 7r°7 with an additional cut of IU < 0.5. a) opening angle 
between the 7r + and the ir°; b) energy of 73. The solid lines are the experimental data 
and the dashed lines are the IB component from the Monte Carlo predictions. 

(±0.1 x 10~ 6 ) was evaluated by shifting the field strength by ±0.2%. The error due to the 
uncertainty in the reconstruction of the kaon stopping x, y position was studied by a shift 
of the kaon stopping point along the charged pion track by one fiber which corresponds 
to 5 mm. The resulting variation of the DE branching ratio was ±0.2 x 10~ 6 . All the 
uncertainties are summarized in Table El 



Table 3: Summary of the systematic errors. 



ERROR SOURCE 


UNCERTAINTY, 
ABr(DE) x 10 6 


7 calibration 


0.5 


7 mispairing 


0.4 


Accidental background 


0.3 


Spectrometer field uncertainty 


0.1 


Kaon stopping x, y coordinates 


0.2 


Total systematic error 


0.7 



5 Conclusion 

We have obtained an improved result for the branching ratio of direct photon emission 
in the K + — > 7r + 7r°7 decay from the analysis of the E470 data. We have extracted 10154 
K + — > 7r + 7T°7 events in the tt + momentum region of 115 to 180 Mev/c, which is about 2.3 
times more than in the previous analysis [16]. In addition to a reduction of the statistical 
error we have also reduced the systematic error by increasing the pairing efficiency of the 
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Figure 7: W spectrum of the observed K + — > 7r + 7T°7 events normalized to the IB compo- 
nent obtained from the Monte Carlo simulation. The solid curve shows the best fit of the 
W spectrum to a sum of the IB and DE components. 

photons. The optimization of the selection cuts allowed us to increase the total number of 
the extracted K + — > 7r + 7r°7 events, and to obtain more sensitivity to the DE component 
part of the K + — > 7r + 7r°7 spectrum. In this analysis we have provided more efficient 
rejection of the K n3 background, which was estimated by MC simulation to be less than 
0.2% as compared to 1.2% in the previous analysis. 

The branching ratio for the DE component in the 7r + kinetic energy region 55 to 
90 MeV has been determined to be 

Br(DE) = [3.8 ± O.B(stat) ± 0.7 (syst)] x 10~ 6 , 

which is consistent with the results of the three last experiments [15-17]. The good 
agreement of this result with the theoretical prediction for the branching ratio of the 
DE component, Br(DE) = 3.5 x 10~ 6 [7,8], supports the hypothesis that the dominant 
contribution to direct photon emission is due to the pure magnetic transition given by 
the reducible anomalous amplitude. It should also be stressed that the result has been 
obtained with the assumption that the INT component is zero. 
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